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ABSTRACT 

Context. Results from gravitational microlensing suggested the existence of a large population of free-floating planetary mass objects. 
The main conclusion from this work was partly based on constraints from a direct imaging survey. This survey determined upper limits 
for the frequency of stars that harbor giant exoplanets at large orbital separations. 

Aims. We want to verify to what extent upper limits from direct imaging do indeed constrain the microlensing results. 
Methods. We examine the current derivation of the upper limits used in the microlensing study and re-analyze the data from the 
corresponding imaging survey. We focus on the mass and semi-major axis ranges that are most relevant in context of the microlensing 
results. We also consider new results from a recent M-dwarf imaging survey as these objects are typically the host stars for planets 
detected by microlensing. 

Results. We find that the upper limits currently applied in context of the microlensing results are probably underestimated. This means 
that a larger fraction of stars than assumed may harbor gas giant planets at larger orbital separations. Also, the way the upper limit is 
currently used to estimate the fraction of free-floating objects is not strictly correct. If the planetary surface density of giant planets 
around M-dwarfs is described as rf/pianct "^ aPda, we find that /3 < 0.5 - 0.6 is consistent with results from different observational 
studies probing semi-major axes between ~0.03 - 30 AU. 

Conclusions. Having a higher upper limit on the fraction of stars that may have gas giant planets at orbital separations probed by the 
microlensing data implies that more of the planets detected in the microlensing study are potentially bound to stars rather than free- 
floating. The current observational data are consistent with a rising planetary surface density for giant exoplanets around M-dwarfs 
out to ~ 30 AU. Future direct imaging surveys will show out to what semi-major axis the above mentioned range of /3 is valid and 
what fraction of the planetary mass objects detected by microlensing are indeed bound. 

Key words. Methods: observational - Methods: statistical - Planets and satellites: detection - Infrared: planetary systems - 
Gravitational lensing: micro - Techniques: high angular resolution 

1. Introduction obj ects are indeed unbound. We summarize the analyses from 

Su mi et all (l201lh in section 2, re-analyze the relevant direct 
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p ■ Gravitational microlensing is a powerful tool to study the initial ^^^^^^^ '^^^^ i„ ^^^ti^^^ 3 ^„^ 4^ ^^^ '^^^^^^^ ^^^ implications 

CN mass function (IMF) of stars and sub-stellar objects for stellar of ^jjg rg-analysis in section 5 

'""I , populations probed along a particular line of sight. It can also ef- 
^ ■ fectively probe the companion mass ratio distribution (CMRD) 

all the way down to terrestrial mass objects. In particular for exo- 2. A large population of unbound planets - 

planet studies the microlensing technique provided some unique derivation Of Current Constraints 
H results in the last years, such as the detection of a cool Su per- 

Earth orbiting an M-dwarf at ~3 AU dBeauheu et al.'. '20061) and Using the datasets from the O GLE and MOA microlensing col- 
unprecedented constraints_on_fliefrec|uencyofcpld Neptune^vlike laborations lSumi et aP (1201 ih identified 10 single microlensing 
an d Jupiter-hke plane ts (ISumi et al.Ll2010tlGould et al.L 120101) . events that lasted only a comparatively short amount of time. As 
I Sumi et aU (1201 ll) reported the detection of an unbound or these short events were not accompanied by a longer signal on 
distantly orbiting planetary mass population based on gravita- which the short event was overlaid the authors concluded that 
tional microlensing results. These pioneering observations could the short events were caused by individual very-low mass ob- 
be interpreted as the low-mass end of the IMF, as bound plan- jects that did not show any indications for a host star. Taking into 
etary mass companions, or as ejected planets. The authors rely account th e detection eflicien cy of different types of microlens- 
on co nstraints from a direct imaging campaign (iLafreniere et al.L ing events ISumi et aP (1201 ih modeled the short events with a 
l2007h to conclude that at least 75% of these objects are not population of Jupiter-mass objects that are LS^^g times more 
bound to any star and were probably ejected from planetary frequent than main-sequence stars with masses <1.0 Mq- The 
systems via dynamical interactions. Here we show that in fact masses of the objects were, however, not directly constrained 
the constraints available from direct imaging are less stringent by the observations and also the existence of host stars could 
than assumed to support the claim that most of the detected only be ruled out to a certain separation. The mean minimum 
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separation between the 10 inferred planets and their potential 
host stars is (amin) ~16 AU but we note that only 2 objects 
have minimum separations exceeding this mean value with ~21 
AU and ~45 AU, respectivel}(j. In order to assess h ow many of 
these p lanetary mass objects might be bound to a star lSumi et alj 
(1201 lb use the results f r om th e Gemini Deep Planet Survey 
(GDPS; iLafreniere et al.L l2007h that aimed at the detection of 
planets through d irect imaging. As explained in their supple- 
mentary material, ISumi et aTl (1201 ll) estimated the upper limit 
for Jupiter-mass planets with separations between 10 - 500 AU 
by taking mean of the 1 Jupiter-mass curve from Fig. 10 of 
ILafreniere eFall (12007), which leads to /„„., < 40% (95% con- 
fidence). As the objects that they detected ha ve minimum sepa- 
rations from their host stars within this range. ISumi et al.l(l201 lb 
conclude that, on average, a star hosts up to 0.4 planets, and 
hence at least 75% of their objects are unbound or free-floating 
planetfl 

3. Re-analysis of Gemini Deep Planet Survey data 

The results presented bv lSumi et al.l(l201 lb imply more than just 
a large population of unbound planetary mass objects. If their in- 
terpretation is correct, i.e., that most of the unbound planets were 
dynamically ejected from their planetary systems, this would 
also require that at least one additional massive planet remained 
in each of those systems implying that, on average, -2.8 massive 
planets formed around each star. As this result would have sig- 
nificant implications on planet formation theory and the overall 
census of exoplanets in the Milky Way, we take a closer look at 
the relevant analysis. We emphasize that we do not re-analyze 
the suggested abundance of planetary mass objects relative to 
main sequence stars, i.e., LS^^g. We rather focus on the con- 
straints that can be derived from direct imaging surveys on the 
10 detected microlensing objects and on the upper limits of stars 
that might harbor planetary companions. 

A very conservative approach would be to read off" directly 
the upper limit for stars having a 1 Mjimiter companion be tween 
10 AU - 500 AU from Fig. 10 ofLa freniereetai] (12007). This 
upper limit is /„„v <92% (with 95% confidence) and is valid for 
any distribution of mass and semi-major axis. 

With access to the data from the GDPS we can directly con- 
strain the occurrence rate of massive planets on orbital separa- 
tions that appear to be the most relevant in context of the mi- 
crolensing results. Using th e data and the approach described 
in detail in ILafreniere et al.l (l2007l) . we computed a maximum 
planet frequency (fmax) pertaining to each of the objects detected 
bylSumi et al. (201 1). We assumed that the planetary masses are 
either exactly 1 Mjupitei or between 0.5 - 3 MjupiteJj and that the 
objects' semi-maj or axes are either e xactly the minimum separa- 
tion flmin given bv lSumi et al.l ( 1201 lb or, to be conservative, twice 
this value. The choice of these two specific separations ensures 
that we have captured extreme cases, i.e., all objects at small or 
all object at large separations. For each of the overall four cases, 
we averaged the upper limits over all 10 objects. The results are 
summarized in Table 1 . 



Table 1. f,„a.x (with 95% confidence) for different combinations 
of assumed planetary masses and semi-major axes. 



' According to lSumi et al.l ( 120111) the minimum separation in AU can 
be computed from the minimum separation in units of the Einstein radii 
of the possible host star rfniin(R£*) via a„,jn « 3AU • c/„,in(R£.)- 

■^ ISumi et sli ( 12011 ) detected 1.8 planets per star from which they 
subtracted 0.4 to correct for those planet planets that could be bound. 
Hence, 1.4/1.8~75% is the fraction of unbound planets in their sample. 

^ This mass range corres ponds roughly to the 68% confidence inter- 
vals that lSumi et al.l i2011l) put on their assumed planetary mass func- 
tion which was centered around ~1 Mjupiter. 
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Jupiter 



0.5 - 3 Mj„pi,e 



78% 
49% 



59% 
29% 



We found that /max can be as high as 78% if we assume that 
the mass of the microlensing planets is 1 Jupiter-mass and that 
the derived minimum separations from the host stars correspond 
to the semi-major axes of the orbits. As expected, the lowest 
value for f,^.^^ is obtained if the mass is assumed to be anything 
between 0.5 - 3 Mjupiter and the semi-major axes are twice the 
minimum separations. We note, however, that this estimate is 
"biased": we implicitly assumed that the planetary mass distribu- 
tion is flat in those cases where we considered the mass interval 
between 0.5 - 3 Mjupiter. However, radial velocity (RV) searches 
for exoplanets suggest that lower mass planets, which are more 
difficult to detect by means of dir ect imaging, are mo re common 
than Jupiter-like analogues (e.g., Mav or et al.l 1201 II) . Thus, this 
upper limit is likely an under-estimate. 



4. First constraints from direct imaging surveys of 
M-dwarfs 

While roughly 3/4 of the stars in the sample from 
ILafreniere et al.l (12007) were FGK stars, it is usually assumed 
that host stars t o planetary objects detected by microlensing are 
M-dwarfs (e.g.. lGould et al.Ll2010b . So, ideally, one wants to do 
a similar statistical analysis based on direct imaging results from 
a large survey of M-dwarf. Recently. iDelorme et al.l (20 1 1 ) pub- 
lished first results from a large deep imaging survey of nearby 
M-dwarfs. They provided the detection limits for 14 objects 
but did not provide a comprehensive statistical analysis of their 
dataset. However, we can use the provided mean detection prob- 
ability of their survey for different planetary masses as a function 
of the separation. We restrict ourselves to separations between 
10 AU - 30 AU, i.e., out to roughly twice the mean minimum 
separation (flmin) derived from the microlensing data. Reading 
off the mean detection probabilities of their survey at 10 AU, 20 
AU and 30 AU for objects with 0.6 Mjupiter and 1.5 Mjupiter and 
averaging them results in an overall mean detection probability 
of (pdetect) ~30%. This is a good approximation for the mean 
detection probability of a 1 Mjupiter object. To be consistent with 
the analysis of the previous section, in particular with Table [Tj 
we also looked at the mean detection probability for the mass 
range between 0.6 - 3.0 Mjupiter for the same separation range, 
which is well approximated by the values for a 1.5 Mjupiter ob- 
ject. Here we find {;?detect) ~45%. 

As no planetary companion was detected in the survey we 
can ask, given these mean detection probabilities, what upper 
limit of stars having a companion in the above-mentioned pa- 
rameter space is consistent with this null-result at the 95% con- 
fidence level. Applying the same Bayesian analysis as done 
above for the GDPS survey and as laid out in Lafreniere et al.l 
(2007) yields upper hmits of /max ^ 59% and /max ^ 40% for 
the 1 Mjupiter case and the 0.6 - 3.0 Mjupiter case, respectively. 
Implicitly, both cases assume a flat distribution of planets in lin- 
ear a space and the second case furthermore assumes that the 
planetary mass distribution is flat as well in the given mass in- 
terval. 
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5. Implications on interpretation of microlensing 
results 

Looking at the upper limits of stars potentially harboring plan- 
ets that we have estimated in the previous sections we find 
tha t in most c a ses th ese values are higher than the 40% used 
by ISumi et aLl ( l2011h . Only in those cases were we looked at 
a planetary mass interval ~0.5 - 3 Mjupiter, with the under- 
lying assumption that the mass distribution is flat, we obtain 
upper limi t s betw een ~30 - 40%. As already mentioned by 
ISumi et al](l2011h . if the objects were Saturn-mass objects rather 
than Jupiter-mass objects, then there would basically be no con- 
straints on their occurrence rate at all. 

It's also worth mentioning that throughout this paper we used 
a confidence level of 95%. Increasing the confidence level would 
also increase the upper limits. 

Finally, it is important to note that the statistical interpre- 
tation of the direct imaging results outlined in Lafreniere et al. 
(l2007h is binomial in the sense that a star either "does not have 
any planet" or "does have at least one planet". Hence, having an 
upper limit /max=50% for stars that may harbor one (or more) 
planets does not mean that, on average, each star has 0.5 plan- 
ets. The upper limits indicate the maximum fraction of stars with 
one or more planets. Thus, regard less of the value ch osen for the 
upper limit, the approach used bv lSumi et al.l (|20I lb to estimate 
the fraction of unbound planets in their sample from the value of 
/inax is not Strictly valid. 

Based on th e points detailed above we believe that 
I Sum! et al] (1201 ll) underestimated the upper limit of stars with 
planets and that, furthermore, the derived upper limit can not 
be used to compute an average number of planets per star. In 
consequence, the main conclusion that 75% of the detected mi- 
crolensing planets are unbound and were probably ejected from 
their planetary systems is overestimated. On the other hand the 
null-result from the planet search around nearby M-dwarfs by 
[Delorme et al. (2011) suggests that it's also unlikely that all of 
the microlensing objects are bound to a star if they are indeed 
Jupiter-mass objects. If there was only one planet of 1 Mjupitei 
between 10-30 AU per star the probability that no object was 
detected around the 14 M-dwarfs is k 6.8 x 10"^. Overall, it 
seems plausible that a certain fraction of the detected planets are 
bound to a star, that another fraction are indeed unbound but rep- 
resent the low-mass end of the IMF, and that a third fraction are 
also unbound but were indeed ejected dynamically from their 
planetary systems. 

6. Predictions for future imaging surveys 

While even more recent d irect imaging survey s for planets 
around solar-type stars (e.g., IChauvin et al.l 1201 Ob were not sig- 
nificantly more sensitiv e to planets i n the relevant parameter 
space than the survey of iLafreniere et al. (2007), future surveys 
on upcoming instruments (e.g., SPHERE, GPI) will shed light 
on the first fraction mentioned ab ove. Also, increasing t he sam- 
ple size of the M-dwarf survey bv lDelorme et al.l (l201lb will be 
of particular interest for the microlensing community. Here we 
can, however, use the upper limits from the current survey and 
combine it with results from microlensing studies to make some 

testable predictions for future surveys^ 

The microlensing results from iGould et all (1201 Ob sug- 
gested that the frequency of planets between 2-8 AU is 
/pianet=0.36±0. 15. In this case, the typical mass of the host stars 
was ~0.5 Mq and the mass regime of the planets was expressed 
in terms of mass ratio q with -4.5 < log{q) < -2 which would 



corres pond to a mass rang e of ~0.02 - 5 Mjupitei. As explained 
above JCassan et all (1201 2b found that the fraction of bound plan- 
ets between 0.5 - 10 AU and with masses between 0.3 - 10 

Mjupiter is /pianet=0.17_Q'gg. 

Using these planet frequencies, we can describe the popula- 
tion of planets by 



dfn-inei — Cirfa^dm da 



(1) 



where C is a normalization constant, m is the mass of the planet, 
a the semi-major axis, and a and /3 are free variables. We can 
now find a value for fi so that the microlensing results are con- 
sistent with the upper limits that we have estimated from the 
M-dwarf survey in section 4. For the power-law slope of th e 
mass function, we use the result from ICumming et al.l (l2008h . 
which in linear mass space gives a - -1.31. This result was de- 
rived from RV searches for massive planets around FGK stars, 
but no equivalent statistical analysis for M-dwarfs is available. 
Integrating equation (1) for varying fi yields different values for 
the constant C, which in turn can then be used to compute the 
planet frequency for the parameter space probed by the direct 



3.0 M 



Jupiter 



and a between 



imaging results, i.e., m between 0.6 ■ 
10-30AU. 

In order to get /pianet 5 0.40 over t his range, we find t hat 
li < 0.61 or;8 < 0.4 9 using the results of ICassan etaP (1201 2b or 
iGould et al] ([2010), respectively. Given the error bars in /pianet 
for those studies, both values for fi are indistinguishable. This 
result shows that current observations do not rule out a rising 
planetary surface density for giant planets around M-dwarfs be- 
tween 10 - 30 AU in linear space. In addition, we note that a 
value for/? in the range of 0.5 - 0.6 provides a direct link between 
the two microlensing st udies mentio ned above and the M-dwarf 
results from [C umming et al.l (l2008b . From RV measurements, 
ICumming et alj(i20 08) estimated fm„„M ~ 0.02 around M-dwarfs 
considering planetary masses between 0.3 - 10 Mjupjter and semi - 
major axes between ~ 0.03 - 2.5 AU (cf. Bo nfils et all 1201 lb . 
Using the approach outlined above and extrapolating these val- 
ues to the parameter space of the microlensing studies yields 
consistent results. Thus, a positive power-law index /3 < 0.5 - 
0.6 is currently able to describe the surface density profile for 
gas giant planets around M-dwarfs between ~0.03 - 30 AU. 

This upper limit for P is, however, in stark contrast to the 
original results from lCumming et aP (l2008b for solar type stars . 
Here, j6 - -0.61 + 0.15 in linear space (see. lHeinze et al.Ll2010b 
for giant planets with periods <2000 days, i.e., a < 3 AU. If we 
take the expectation valu e 13 - -0.61 and extrap olate again the 



M-dwarf fr equency frorn Curnmin g et al.l (120081 tojhe parame- 
ter space of lGouldet"an (l2010b and Cassan et aL| (1201 2b we find 
that /planet ~ 0.05 and /pianet ~ 0.03, respectively, which is not 
consistent with observational results within the given error bars. 
Taking into account the 3-cr error bar and using /? - -0. 1 6 would 
result in /piane t ~ 0.11 a nd /pianet ~ 0. 06, again out side the error 
bars given by Gould et al. (2010) and lCassan et al.l (12012). 

Future surveys with SPHERE and GPI, or an extension of 
the survey done bv Delorme et al. ( 2011) . will reveal the number 
of bound giant planets between 10-30 AU around M-dwarfs, 
empirically constrain the power-law slope /3, and test to what 
separation the predicted value range is valid. These surveys will 
provide some further hints as to what mix of bound and unbound 
planets is the most likely - and les s extreme - explan ation for the 
microlensing results presented bv lSumi et al.l (1201 lb . 
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